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Abstract 

We studied the microscopic mechanism of multiferroics, in particular with 
the "spin current" model (H. Katsura, N. Nagaosa and A. V. Balatsky, Phys. 
Rev. Lett. 95, 057205 (2005)). Starting from a system with helical spin con- 
figuration, we solved for the forms of the electron wave functions and analyzed 
their characteristics. The relation between ferroelectricity and Dzyaloshinskii- 
Moriya interaction (I. Dzyaloshinskii, J. Phys. Chem. Solids 4, 241 (1958) 
and T. Moriya, Phys. Rev. 120, 91 (I960)) is clearly established. There is 
also a simple relation between the electric polarization and the wave vector of 
magnetic orders. Finally we show that the bond-bending existing in transition 
metal oxides can enhance ferroelectricity. 
PACS: 75.80.+q, 77.80.-e 
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0.0.1 1. Introduction 

Experimental findings [1-4] in recent years have revived the interest in multifer- 
roics. They showed that magnetic and ferroelectric orders are closely related[5- 
8] . What is more intriguing is that only certain types of magnetic orders, namely 
helical spins and frustrated spins, can be coupled to ferroelectricity[9]. It is 
this fascinating interplay between ferroelectric and magnetic orders that has 
attracted many researchers. There were already models based on the Ginzburg- 
Landau theory[10,ll] that provide instructive physical description of the sys- 
tems. As for the microscopic mechanism, there are currently two schools of 
theories. One of them proposed that the electric polarization and the anomaly 
of dielectric constant come from atomic displacements. The displacements or 
phonons are in turn, coupled to spins[12-14]. Though proposed for systems of 
orthorhombic structure, it is more readily applied to multifcrroics of hexago- 
nal structures, such as HoMn03, as there is experimental evidence of atomic 
displacements from neutron scattering data[15]. The second school of theory 
proposed a new possibility: electric polarization coming from electronic wave 
function and thus density distribution. Katsuraet. al.(KNB)[16] predicted that 
the magnetoelectric effect can be induced by "spin current" [17]. The coupling 
between "spin current" and internal electric field has the same form as that of 
Dzyaloshinskii-Moriya interaction (DM) [18,19] or AC-effect[20] where the mo- 
tion of a magnetic moment is coupled to electric field. In this latter theory, the 
atomic displacement is not essential. On the other hand, spin-orbit interaction 
is indispensable in generating electric dipole moments. 

The " spin current" model, though a bright idea, needs additional substanti- 
ation in order to be applied to physical systems. Jia[21] et. al. gave a detailed 
calculation of this model. Their results showed that the "spin current" model 
is able to explain at least semi-quantitatively many experimental data. This 
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model was also applied to the systems with e g orbitals such as TbMnC>3 [22] . In 
our opinion, the foremost task is the embodiment of this idea in a crystal in 
which completely different properties and behaviors can emerge from complex- 
ities and interrelations between electrons, spins, and lattice structure. Equally 
important is the calculation of the magnitude of electric polarization induced 
by "spin current". According to KNB, the polarization is of the orders el(t/ A) 
or eI(t/A) 3 for one or two holes. Here / is the expectation value of length, t 
the hybridization energy and A the energy difference between d-orbitals and 
p-orbitals. According to their estimation, t = V{pdit) <~ O.leF, A 2eV, 
and el /a 3 ~ 10 4 fiC/m 2 where a is the lattice constant. One can see that the 
magnitude is compatible with experimental data only in very favorable condi- 
tions. Therefore, it is desirable to conceive possible and realistic mechanism to 
enhance the ferroelectricity-magnetism coupling. 

KNB derived an elegant expression for the electric polarization of a three- 
atom system P <~ e/ei2 x (ei x e<i) where ei and e*2 are the directions of spins 
of the transition metal ions and e\i is the bond direction. We shall see how 
the expression conforms in a crystal and how it is related to the wave vector of 
the helical magnetic order. We will explain the reason why it is advantageous 
to have the helical spin configuration for ferroelectricity and also elaborate the 
roles played by the spin-orbit interaction. In fact, we will show that the electric 
polarization comes directly from DM interaction. Finally, we point out that 
a common feature in transition metal oxides, the bond-bending, can enhance 
electric polarization. 

0.0.2 2. Description of the system 

We are going to consider two features in the system, the helical spin configura- 
tion and bond-bending. Mostovoy[23] studied a system of degenerate double- 
exchange interaction and next-nearest neighbor hopping. He found that un- 
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der certain conditions the helical spin configuration is stable. Thus, we shall 
take that as our starting point. The transition metal ion, with position vector 
Rj.m = Rj+ ~r *m where Rj is the position vector of the j-th lattice point and 
~r m is the position vector of the m-th ion in the basis, has the following form 
for its spin (presumably those of t^g electrons) 

S j, m = S[(e x cos 4> + e y sin </>) sin(~<f • R + e z cos(~q ■ R jiTn )] (1) 

The helical spin order has a wave vector ~q. As shown in Fig. 1, the projection 
of spins on xy-plane makes a fixed angle </> with x-axis. 

Hund's coupling in the transition metal ions, — Jjy£ S j,m • ~Sj.m, is the dom- 
inant mechanism in the system. If the local spins in cq. (1) are treated as 
classical spins, their effect on the spins of hybridizing electrons, denoted by 
~~Sj, m , is equivalent to an effective magnetic field. Thus, the cigen state of the 
spin of a hybridizing electron on the m-th transition metal ion with lower energy 
is 



where 0j : m = ~Q • Rj,m/2- This will affect the hybridization of orbitals. 

Next, we consider the hybridization. KNB studies a TM-O-TM three-atom 
triad. Due to the symmetries of orbitals, the hybridization arises from ir- 
bonding. If the TM-O-TM triad is not linear (illustrated in Fig. 2a where 
a is the bond angle) , the p x orbitals of oxygen (if the bond is approximately in 
the x-direction) can also take part in electron transfer, cr-bonding can be realized 
and the hybridization energy can be greater. In fact, the bond-bending occurs 
quite often in transition metal oxides. See, for example, reference 24. The ab 
plane of the crystal is shown in Fig. 2b where the solid dots and circles denote 
the transition metal ions and oxygen atoms respectively. There are two distinct 




(2) 
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transition metal ions in the basis, and thus m = 1,2, due to bond-bending. 
As a result, two transition metal ions and four oxygen atoms form the basis of 
the crystal which is enlarged by bond-bending. The atoms are labeled so as 
to facilitate later deduction. As mentioned above, bond-bending will affect the 
hybridization energies. They can be found in Slater and Koster[25]. 

Bond-bending has other profound effects. It changes the symmetry of the 
surrounding of atoms and thus the symmetry of crystal field. In an orthorhombic 
crystal (a ^ b) with bond-bending, the xy-orbitals of the transition metal ion 
will mix with the x 2 — y 2 orbitals. The resulting orbital has the form 

cos fi\xy) Tsin/3|a; 2 -y 2 ), (3) 

where the ± sign is determined by the direction of the displacement of the 
oxygen atom away from the line joining two transition metal ions. (3 and a are 
of the same order of magnitudes. Ideally, if the crystal field is determined by the 
four nearest oxygen ions on the ab plane, then [3 = a' + 5a' where a' = it — a. 
The deviation 5a' is due to the fact that the symmetry of crystal field is not 
determined by bond-bending alone. We expect 5a' to be finite but smaller than 
a in orthorhombic crystals of manganites but large in compounds like NI3V2O8. 
However, to get its magnitude, one need to carry out a very precise first-principle 
calculation which is beyond the scope of the current paper. 

The hybridization of the p-orbital of oxygen atoms and the xy-orbitals of 
transition ions is equal to ±(v / 3/2) sin a sin(a / 2)V (pda) where the sign is again 
determined by the direction of the displacements of oxygen atoms. That of 
the p-orbital of oxygen atoms and the x 2 -y 2 -orbitals of transition ions is equal 
to — (v / 3/2) cos a sin(a/2)V (pda). Although the x 2 -y 2 -orbitals have greater hy- 
bridization energy, its amplitude is smaller as one can see from (3). The re- 
sulting hybridization energy is ±(\/3/2) sin(o: / — /3) sin(a / 2)V (pda) . We note 
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in passing that since a is close to n, if \(3 — a'\ w 7r/12, then (\/3/2) sin(a' — 
/?) sin(a/2)V(pda) « 0.2F(pd(r), which is often greater than T^(pd7r). 

0.0.3 3. Calculation 

Due to the Hund's coupling, the spins of the hybridizing electron are projected 
to the local spins and the hybridization energy is modified. It is the main part 
of our Hamiltonian. We considered a bond only. Hence, the orbitals involved 
are that in (3) of transition metal ions, p x — orbitals of oxygen atom 3 and 4 and 
p y — orbitals of oxygen atom 5 and 6 where the numerals are shown in Fig. 2b. 
Later, we take the spin-orbit interaction, and more orbitals, into account. It will 
be treated as a perturbation. Our Hamiltonian thus has two parts, H — H +Hi : 



n.n. 

(4a) 

and 

Hi = A^ I jtTn ■ ~Sj, m (4b) 

j,m 

where the hybridizations energy V = (\/3/2) sm((3-a') sin(a/2)V(pda), is only 
among nearest neighbors. A is the strength of spin-orbit coupling. e p , £d and 
c pil, t(J.)i c djm are the energies and field operators of the p— orbitals in oxygen 
atoms and d— orbitals in transition metal ions. / and m are the indices of oxygen 
atoms and transition metal ions in the basis. As shown in Fig. 3, m = 1, 2 for 
manganese ions and I = 3 — 6 for oxygen atoms. As previously mentioned, the 
localized spins were approximated as an effective magnetic field of angles 6j im 



6 



and <f>. We now made the following transformation 

p a ^ = cos(~q ■ Ri,i/2)cpU,-\ + e 1 ^ sm(~q ■ Ri,i/2)cpH,i (5a) 
Pa, I = cos(~q ■ ~Rij/2)c pi i,i - e~ %4 ' sin(lf • ~R lt i/2)c p a A (5b) 

where Rij = Ri+ ~r^i with ~ri being the position vector of the i-th oxygen 
atom in the basis without bond-bending. The Hamiltonian becomes 

n.n. 

- Vl (-l) m e^ sin Se^pu^+H.c. (6) 

where d jm = c djm , S6 l = a q x /2, rji = 1 for I = 3, 4 and 56i = b q y /2, r]i = -1 
for I = 5,6. a and bo are, respectively, the distance between transition metal 
ions along x and y direction. The Hamiltonian in momentum space is 

Ho = Y, E pP\k,*P*,° + E ^ d lkdmk £ V[(-l) 1 cos 69^ ■^'■•-^■™>4, fc p Ifc , T 

n.n. 

-mi-ire^sinSete^^-K^Ul^k^+H.c.. (7) 
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The following is its matrix form: 
/ 



V 



£d 





V +X c 3 


-V +X e" 


b s 3 


-V- x c A 


-V- x e^s A 


V-yC 5 




-V+yCQ 





£d 


V- X C 3 


V- x e +i < 




~V +X C4 


V +x e^ S4 , 


V+yC 5 


-V+ye^SB 


-V- y c 6 


V- X c 3 


V +X C 3 


£p 






















-V- x e-^8 3 


V +X e~^s 3 





Sp 





















-V- X c 4 










£p 














-V +x e-^ Si 















Sp 











V+yC 5 


V- y c 5 
















Sp 








V +y e~^s 5 





















Ep 





-V-yC 6 


-V+yCe 






















Sp 




-V+ye-^se 


























(8) 



Here the order of matrix elements is m = 1,2, and then I = 3 |,3 J., 4 j,4 j 
etc.., V± x = Ve ±lk * aa , V± y = Ve ±lk y b ° and = sin 89 l , c, = cos^. It can be 
diagonalizcd and the eigen values are 



£d - £ P 



f+W 2 ±V 2 \Y^ cos 6e nn e 2i ^-^r 



(9a) 



E- 



2± 



r £d — £i 



) 2 + w 2 ±v 2 \Y^ 



cos56 nn e 2ik - R r>n\ (9b) 



E 3 



(9c) 



where e p is six- fold degenerate and R nn and <!>#„„ = 258i are the position vector 
and spin angle difference between two nearest neighboring transition element 
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ions. The eigen vectors are 



cos 8 



(10a) 



e' v sin( 



(10b) 



^(-l)^' 1 ^ 

E2+ — £ P 



.V(-l)'e 



i k ■ r 



E2- — £ P 



-il> v i{~r-Rj,i)\( 

(10c) 
-ip p i( 1 r-'Rj,i) ( 



cos 0,- 



e' v sin 1 



cos( 



e v sintf,- 



(lOd) 



where A m u2\± are normalization constants, ~r ; m = V; — "r m , and 



A m l± _|_ n il. 

A 



E cos 86ti 



2ik R r 



m2± 



E cos 50 n 



o2ik R r 



(11) 



Note that ip p i = ip px for I = 3, 4 and = for I = 5,6. The eigen vectors 
of e p do not concern us because they are non-bonding states. If there is no 
bond-bending, then there are only 7r-bonding for xy-orbitals. It turns out that 
the eigen vectors have very similar forms as those in eqs. (9). They are shown 
in Appendix. 

0.0.4 4. Spin-orbit interaction and polarization 

We now introduce the spin-orbit interaction in eq. (4). Its effect can be ex- 
pressed in the following relations 



I ■ s \xy, t) = i\zx, I) + \yz, [) - i\x -y , |) 



(12a) 
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I ■ s \xy, 4) = i\zx, t) - \yz, 1}+i\x - y , |) 



(12b) 



~s \x 2 - y 2 , T> = \zx, I) - i\yz, J.) + i\xy, T) 



I -~s\x 2 - y 2 ,[) = -\zx,1) - i\yz,1) - i\xy,[). 



(12c) 



(12d) 



Therefore, if the spin-orbit interaction is treated perturbatively, the wave func- 
tions tpd(~r — Rj. m ){cos8j. m \ |) + e 1 ^ s\nOj. m \ [)) in cqs. (10) will be replaced 

by 



cos 



COS 8j. r , 



sin V jtT , 



AE, 



cf 



/ -icoaO jtm x 

Wx 2 -y 2 { r - Rj,m)[ J 

ie 1 ^ sin 9i m 



oi<t> 



_> — > / te 'sin^ m \ _> — > / — e^sin^j iTO \ 

+ ^x(r - Rj, m )[ J +1pyz(r - Rj,m)[ )}} 

icos8j m 



cost 



3,m 



r y — » / COS "j rn \ A f y — > ( 

e 1 ^ sin 0.- ™ c ^ 



i cos # j >m 



-ie 1 * sin 4 



COS( 



— ie v sm( 



-i cost 



(13) 



where AE c f is the energy difference between cos/3|xy) =F sin/3|x 2 — y 2 ) in (2) 
and sin/3|xy) ± cos/3|x 2 — y 2 ). Substituting (13) into eqs. (10), we found that 
the polarization per unit cell is 



k J 

8V2sm/3VX 



7r(e P - e d )AE cf 



(epe z )[cos(j)sm(a q x /2) + sin0sin(6 o %/2)] (14) 
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where 



P « - j d 3 !^ ip* zx {~r)ziP px {~r - a e x /2) w - ^ d 3 ^* z ("r )x^ p2/ (T > - 6 e y /2) 

(15) 

with a and 60 being the lattice constants without bond-bending. Another way 
of writing them should be "a = (Rj+i ~ an d ^ = {Rj+i — -Rj)/2 for 

site j + 1 situated at x- or y-direction away from site j. There are ten bands in 
our calculation. The main contribution comes from the topmost occupied band. 
We have also made the approximation E n+ w e p , taking the advantage of the 
fact that V <C \s p — ea\. If even number of bands are occupied, the polarization 
will be much smaller. The polarization produced by two bands tends to cancel 
each other. As a result, there is an extra factor of 4V 2 /(e p — Ed) 2 which comes 
from the denominator of eqs. (10): 



32V2sin/?y 3 A 
7r(e p - e d ) 3 AE cf ' 



P-$ ~ — „. V 3af (epe z )[cos0sin(Q O (7 z /2) + sin<^sin(6 %/2)] (16) 



To see more clearly how polarization and q are related, let Up be the space 
inversion operator. We then have 

^,^) = v>_ t ,-^)- ( 17 ) 

Under the space inversion, the displacements of oxygen atoms and hence, a, [3 
and V change sign under inversion 

P ? = -e (UpUpVUpUp) = -P_^ (18) 

and the polarization is an odd function of ~~q . The form of sine function in eqs. 
(14) and (16) seems to be a natural form. In the continuum limit (oq w b — > 0) 
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or the long wavelength limit (~q — ► 0), we found 

P-j ~ 0.01epaoe z (q x cos <j) + q y sin (f>) ~ 0.01epao~q x /i (19) 

where h is the helix axis unit vector. We have assumed that q z = since our 
calculation was performed on xy-plane. If one recognizes that the true meaning 
of ei2 is the spatial direction alone which the spins propagate, i.e., ~q, he can 
see how the form P ~ e/ei2 x (ei x e 2 ) given by KNB can be transformed in 
the presence of magnetic orders. 

0.0.5 5. Discussion 

In order to see how magnetic orders and electric order are coupled, let us go 
back to the original work of Moriya[19]. He derived the following expression: 

Hdm = n,n-i ■ (S n x S n-i) (20) 



N 



where 

_^ J(n, n', m, n') ( n\ I j|m ) J(n, n',n, m') (n'\ I j\m'\ 

d n i = i\Y ^ '- - i\Y ^ '-. 

(21) 

Here, J(n,n' ,m,m') is the exchange interaction strength and / , denotes the 
angular momentum of the electron at site i. 

In our starting Hamiltonian, the exchange interaction comes from the charge 
transfer energy e p — Sd and hybridization energy V[27]. Combined with the spin- 
orbit interaction in eq. (4b), DM interaction is clearly present in the system we 
considered. We can recast the wave functions we got previously in a form similar 
to Moriya's by treating the hybridization energy and spin-orbit interaction as 
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perturbations, 



_> (m\H v \n)X(n\l\l)-(a n \T\ai) 

P* = -e£ (e d -e p )AE cf {1,<ti\ r \m,a m ) + ex. (22) 



where H v is the hybridization energy (the third and fourth terms of eq. (4a)). 
\l) and | ct;) are respectively the spatial and spin part of an intermediate state. 
|n) and |^) arc states of d-orbitals and \m) p-orbitals. (n\Hv\m) /{ed — e p ) is 
the exchange part with Hv also containing information of spins. A (l\ I 
(ai \~s\<r n ) is the spin-orbit coupling. Eq. (22) is also applicable to the situations 
without bond-bending. In that case, the 7r-bonding will be considered as they 
were in KNB's original work and Hy is the third term of eq. (A-l). One can see 
that eq. (22) has the same origin as Hdm in eqs. (20) and (21). However, we 
have to note that the helical spin configuration is not caused by DM interaction 
whose strength is too small. Rather, it can be due to the next-near-neighbor 
hybridization as shown in ref. 23. 

It is easier to analyze with the following form: 

p (m\H v \n) X(n\e hij r h pi\l) {a n \sj\ai) (l,<Ti\r k \m,a m ) 

? < fc ~ 6 ^. (e d -e p )AE cf ^ 

(23) 

where Ehij is the antisymmetric Levi-Cevita symbol and r k is the fc-th compo- 
nent of ~r in space. We consider the mirror symmetry of above equation. Hy 
may change sign under mirror reflection operation because of the orbital wave 
functions involved. It is also related to the direction of the displacements of 
oxygen atoms. For example, the (pd7r) part of E y _ xy and (pder) of E x x 2_ y 2 ( in 
Slater-Koster notation) with bond along x-direction change sign if one makes 
the operation x — ► —x. On the other hand, the (pder) part of E x<xy (due to 
bond-bending) does not change sign if one makes the operation x — ► — x or 
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y — ► —y. In the previous section, we calculated the polarization of a planar 
crystal. Now we consider a more general case of orthorhombic structure and 
helical spin configuration. It is easier to catch the essence if one considers the 
(pd7r) part of E ViXy or the (pda) part of E XjX 2_ y 2 of Hy, which is applicable to 
KNB's original work and to the x 2 — y 2 orbital part of our work respectively. 
Assuming the bond direction of Hy is in the m-direction then polarization is 
finite for either k = h, m = i or k = i, m = h. As a result, eq. (23) can be 
simplified as 



where Hy,i denotes the hybridization bond along i— direction and i, j and k are 
cyclic. 

The spin part needs more attention. (<j n \sj\ai) give rise to a spin state 
\<7i) different from the original state \a n ) and (<Ji\<J m ) gives rise to a interesting 
contribution. Only the imaginary part needs to be considered because (Z|Zj|n) is 
imaginary: 

lm({a n \s x \ai) (<ri, x \a n )) = sm4>sm(~q • i?jv,JV-i/2) = (~s n x ~~Sn-i)\x (25a) 
lm({a n \s v \(7i) (cri,y\(Tn)) = - cos^sin("g • R n ,n-i/2) = (~sn x ~s N -i)\ y 



for the spin configuration in eq. (1). Here the additional subscript j of \<jij) 
denotes that it comes from (a n \sj\ai). The right hand sides of eqs. (24) are very 
similar to respective components of S n x S n-i which appears in Hdm- The 
difference is in the arguments of sin functions. The factor 1/2 arises because in 
our model it is the hybridization electrons that mediate the exchange interaction 
while in the original DM interaction it is direct exchange. Substituting eqs. (25) 




(m\H Vi i\n) \(n\lj\l) (a n \sj\ai) (l,ai\r k \m,cr m ) 



(24) 



(25b) 
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into eq. (24), we found that 

Pf, k = -2e £ ^^(n|y/)(^x^_0|,^ + c.c.. (26) 
, . (£d — £p) £±&cf 

The connection between the electric polarization and spin current can now 
be established. A common definition of the spin current is 

A = i^X^Riv+aMRiv - ^O-jdRw+aJ^li ( 27 ) 

where 1? is the lattice vector of length a in the direction of i and in our case 
t w V 2 /\e<ip\- Note that the electrons hop along the direction of ~a. Equation 
(27) can be derived by discretizing the form ih[i)) + ajdi^ — (diil> + )ajtp\/Am. The 
form in cq. (27) manifests itself in helical spin configuration. We calculated its 
expectation value with the state in (2) and found that 

3x) = ^cos0sin(g x a o /2) = - ^C«iv x ~t N -i)\ y , (28a) 

iy) = -^sm<j>sm(q y a /2) = -^jrC?iv x ~s N -i)\x- (28b) 

The results are same as those in eqs. (25). Introducing the SU(2) vector poten- 
tial^] A l j = —eijkhEk/Am e c 2 where m e is the electron mass, the DM interac- 
tion can be rewritten as 

H DM =eJ2A]jf. (29) 

where 

A ^ = -^AE7 f eB ^ (30) 

is the guage field coupled to the spin current. It comes from the term E ■ (~p x ~a) 
in Pauli's equation. See for example, ref. 28. The polarization results from the 
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perturbation of Hdm'- 

E <M|Erf> (0WM) (31) 

M E « M 

where |0), the ground state, has a component (V/ed P )\tp P )- If we take E m to 
be the exchange energy J, then eq. (25) and eq. (30) are equivalent in view of 
cqs. (26) and (28-30). 

Now we can see more clearly what the origin of internal electric field Ek is. 
From eq. (30), we found 

— > 8m e c 2 JX ^ — > ,„„. 
E = 7 a j-i e B x I 32 
eta AE c f 

where the factor m e c 2 will be cancelled by its inverse in A. The small factor J/t 
arises because of the cancellation between different band. See also eqs. (14) and 
(16). The electric field originates from the perturbation of spin-orbit interaction. 
The spin-orbit interaction changes the angular dependence of the d-orbital wave 
functions and hence, the electron density, via the hybridization of d-and p- 
orbitals. This effect can be interpreted as the result of an internal electric field. 
However, the electric field thus created does not necessarily give rise to net 
electric polarization. Certain environments are more advantageous than others. 
As one can see from eqs. (20), (25) and (28), the helical spin configuration is 
apt to provide spin current, net electric polarization and multiferroics. 

The magnitude of polarization is also important. If there are odd number 
of filled bands then the polarization is of the order P-^/fl where f2 <~ 250 A 3 
is the volume of a unit cell and eq. (14) is used. If we take \e p — Ed\ ~ 2eV, 
AE cf ~ 2.0eV, V = (V3/2)sin(a - 0) sin{a/2)V(pdcr) - 0.2eV, A - 0.05eV 
and sin/3 ~ 0.1 then P <~ 10yuC/m 2 for p ~ 0.1A The bond-bending activated 
hybridization V is in general greater than V{pdir). and psin0/I where / is 
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defined in eq. (A-6), is of the order V(pda) sin ^/V(pdn) <~ 1 for /? 7r/12. 
Hence bond-bending gives larger polarization in many oxides. For example, 
a ss 5-7I-/6 in manganites[26] and in compounds such as Ni 3 V20s the bond 
angle of Ni-O-Ni a can be as small as 7r/2[6]. Hence, a and f3 may both be 
large and the environment is favorable to ferroelectricity. Furthermore, p <~ 
0.1A can very well be an underestimation because the 2 ~ has a much larger 
radius than a neutral oxygen atom. All things considered, P can be an order 
of magnitude greater than previously estimated. Atomic displacement can also 
enhance electric polarization by destroy cancellation. However, if it has its own 
wave vector and it is not commensurate with ~q then there is no net polarization. 

In conclusion, we have analyzed the conditions for the emergence of fer- 
roelectricity due to magnetic orders. We found a simple relation between its 
wave vector ~q and polarization. Furthermore, the physical picture of the cou- 
pling between magnetic orders and ferroelectricity is made clear. Multiferroics 
is created by a generalized version of Dzyaloshinskii-Moriya interaction in the 
environment of certain spin orders, preferably helical. Above findings can also 
be applied to systems without bond-bending but with 7r-bond hybridization. 
The bond-bending tends to enhance polarization and it may be important for 
certain compounds. 

The author benefited from the activities of "quantum novel phenomena 
in condensed matter" focus groups of NCTS, Taiwan and discussion with S. 
Maekawa. This work is supported in part by the National Science Council 
under the contract NSC 95-2112-M-002-048-MY3. 

0.0.6 Appendix 

In this Appendix, we present the eigen values and eigen vectors of the 7r-bonding 
systems. It is simpler because there are only one transition element ion (thus 
the dropping of the index m) and two oxygen atoms in the basis. We consider 
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a planar crystal on the xy-plane. For bonds along x(y)-direction, p y ( x )-orbital 
of the oxygen atoms and the xy-orbital of the transition element ions form ir- 
bond. The p z -orbitals are ignored because they do not give rise to polarization 
in z-direction.orbitals are considered. Thus, we have the Hamiltonian: 

H o = J2 £ P c U,<y c P ll ^ + J2 £ d c dj c dj-J2 V'(-l) l [coa Ojc^Cpii^+e^ sin 6jC ] dj c p i hl ]+H.c 

n.n. 

(A-l) 

where V = V(pd,Tr) and I = 1,2 for the oxygen atoms on x-axis and y-axis 
respectively. Making a transformation similar to that in cqs (5), we obtain the 
Hamiltonian in momentum space: 

H = $> P j& i<r Pifc, ff + £> < 4d fc - ^'(-l^cos^e^- ^dj^M 

n.n. 

-e^smSO^^-K^dlp^l+H.c.. (A-2) 
Now 661=3,4 = q x ao/2 and SOi—^q — q y bQ. We can solve for eigen values 



E[ ± = £ A±£p + /(£^iP)2 +4y /2 ±y /2|^ cos2 ^ e 2ife.^ 1 | ( A _ 3a ) 



E' 2 ± = £1±£p - ( £ -±- ^)2+4T/' 2 ±y' 2 |^cos2^e 2 ^-^i| (A-3b) 



E — e r . 



(A-3c) 



where e p is two-fold degenerate. The eigen vectors are 



. E' 1+ s p 



cos ( 



e v sm 



(A-4a) 
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fe • r i 



cost 



e** sin ^ 



(A-4b) 



^2+ C P 



i>pi(~r-Rj,i) ( 



cost 



sin 0,- 



(A-4c) 



j i 2 ~ p 

(A-4d) 

where A' 1± and A' 2± are normalization constants. One can easily see the similar- 
ity between the wave functions with or without bond-bending. The polarization 
can be computed as eq. (13). The result is 



cos 8j 
sin 6j 



?u = <_ e r> « ^ 



— (eJe z ) [cos sin(a g , x /2) + sin sin(6 %/2)] 

(A-5) 



where Ai?^ is the energy difference between t2 g states and e g states and 



/ = d d r iP* yz ( r )z^ py { r -a e x /2) = d 6 r i>* zx { r )zip px ( r -b e y /2) (A-6) 
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Figure caption 



Fig. 1 Schematic helical spin configuration. cj> is the angle between the 
projections of spins on xy-planc and x-axis. The angle between spins and z-axis 
is twice of 9j !m = ~q ■ Rj iTn /2 defined in text. 

Fig. 2 (a) Bond angle a. (b) A two-dimensional lattice with bond-bending. 
The solid dots and circles denote the transition metal ions and oxygen atoms 
respectively. The atoms are labeled so as to facilitate later deduction. 

Fig. 3 Polarization P z versus (j> and <p q where 4> q is the angle between ~~q and 
x-axis. 

Fig. 4 A spin current along the Mn-O-Mn bond (y-componcnt of spins 
moving in x-direction.) It can be coupled to the z-component of electric field 
and induce polarization P z . 
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Fig. 2a 
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Fig. 2b 



25 




2G 



Fig. 4 
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